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SECTION  I 
INTRODUCTION 


Present  implementation  of  the  beamformer  for  the  BQR-7  submarine  sonar  provides 
for  stave-to-stave  delaj'S  only  .'^Consequently,  it  lacks  the  capability  of  steering  or  broad- 
ening the  vertical  beam  by  utilization  of  a time  delay  taper  across  the  vertical  dimension  of 
the  array.  In  esseiicja^'-i^ie  array  is  optimized  to  receive  target  radiated  noise  through  the 
surface^dupt'inbde.  For  short  ranges  between  the  target  and  the  receiver,  this  gives  satis- 
factory operation  since  most  of  the  energy  is  in  the  surface  duct. 

longer  ranges,  however,  most  of  the  energy  may  be  propagated  via  the  bottom  bounce 
mode.  This  results  in  less  than  theoretically  maximum  performance  for  long  ranges  since 
the  energy  enters  through  the  sidelobes  of  the  vertical  beam  pattern.  For  targets  of  interest, 
the  signal-to-noise  ratio  (SNR)  at  short  ranges  is  high  enough  so  that  one  might  consider 
sacrificing  some  system  gain  at  short  ranges  if  the  detection  performance  at  longer  ranges 
can  be  significantly  improved. 

There  are  several  methods  which  might  be  utilized  In  upgrading  the  performance  at  long 
ranges.  One  method  consists  of  steering  the  beam  downward  to  intercept  the  energy  arriving 
via  the  bottom  bounce  mode  through  the  main  beam  of  the  vertical  pattern.  An  alternate 
method  consists  of  broadening  the  vertical  beam  pattern,  without  vertical  steering,  to  inter- 
cept both  bottom  bounce  and  surface  duct  energy  with  the  main  beam.  A third  choice  consists 
of  a suitable  combination  of  vertical  steering  and  broadening  to  optimize  the  detection  per- 
formance for  all  ranges. 

A study  was  undertaken  tovmvestigate  the  various  aspects  of  vertical  beam  steering  and 
broadening  and  its  effect  on  detecthm  performance.  An  attempt  was  then  made  to  verify  the 
theoretical  results  by  beamforming  actuql  sea  test  data  previously  recorded  at  the  hydrophone 
level  aboard  an  SSB(N)  640  class  submarin^s^ 

The  feasibility  of  designing  a compromise  beamformer  to  handle  numerous  classes  of 
submarines  was  also  examined.  Specifically,  the  classes  of  submarines  examined  were  the 
SS(N)  594,  SSB(N)  598,  SSB(N)  608,  and  the  SS(N)  671. 

The  results  of  the  vertical  beam  coverage  study  as  well  as  the  compromise  beamformer 
feasibility  study  are  reported  herein. 
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SECTION  II 
ASSUMPTIONS 

1.  DETECTION  MODEL 

The  analyi.ical  model  used  to  investigate  the  detection  performance  of  the  SQR-7  DIMUS 
beamformer  calculates  the  SNR  at  the  output  of  the  detector.  In  general,  the  SNR  at  the 
output  of  a passive  detection  system  is  given  by 


SNR  = 


where 


P (f,  R)  = signal  power  spectrum  at  the  input  to  the  detector  as  a function  of 
frequency  (f)  and  range  (R) 

P^(f)  = noise  power  spectrum  at  the  input  to  the  detector 

T = the  integration  time  of  the  output  filter 
fj^  = lower  frequency  limit  of  receiver  operating  band 

f_  = upper  frequency  limit  of  receiver  operating  band 

The  signal  power  spectrum  at  the  input  to  the  detector  is  calculated  from  the  following 


expression: 


^ L^(f)  R)  Y(f) 
,(t,  R)  Nw,(f.  R) 


where 


L.p(f)  = target  radiated  noise  power 

Np  = number  of  paths  from  the  target  to  the  receiver 

Ng.  (f,  R)  = vertical  beam  shape  factor  (the  loss  suffered  by  the  signal  arriving  from 
^ the  i^h  path  due  to  arrival  off  the  main  response  axis  of  the  array  or  due 

to  beam  broadening) 

Tp(f)  = the  transducer  power  response  characteristic 
Y(f)  = the  pre-emphasis  filter  characteristic 
N\y.(f,  R)  = the  propagation  loss  for  the  i ” path  from  the  target  to  the  receiver 
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The  noise  power  spectrum  at  the  input  to  the  detector,  on  the  other  hand,  is  given  by; 


where 


P (f) 
n'  ’ 


* '^ON«]  ■'■r® 


(3) 


Loxt(^)  = own  ship's  self-noise  level 
Loj^(f)  - ocean  background  level 
Npj(f)  = directivity  index 


A computer  program  has  been  developed  to  calculate  the  SNR  as  given  in  Equation  (1) . 
The  program  calculates  the  signal  power  spectrum  and  the  noise  power  spectrum  at  discrete 
points  within  the  receiver  operating  b|.nd  by  use  of  Equations  (2)  and  (3).  Equation  (1)  is 
then  evaluated  by  numerical  integration.  The  computations  are  repeated  for  each  of  the 
desired  ranges. 


Since  the  objective  of  the  present  study  was  to  determine  the  effect  of  vertical  beam 
broadening  and  depression,  the  previous  computations  were  carried  out  under  numerous 
combinations  of  vertical  depression  angle  and  vertical  beam  broadening.  The  effects  of 
broadening  can  be  accounted  for  by  appropriate  adjustment  of  the  vertical  beam  shape  factor. 

The  directivity  index  for  the  BQR-7  DIMUS  utilized  for  the  computations  was  taken  from 
a rough  draft  copy  of  a memo  (no  title)  issued  by  USNUSL,  dated  31  July  1967,  No  change  in 
directivity  was  assumed  due  to  beam  depression.  However,  adjustments  were  made  in  the 
vertical  beam  shape  factor  to  account  for  the  loss  in  directivity  due  to  beam  broadening. 


The  vertical  beam  shape  factor  was  assumed  to  be  that  of  a three-element  line  array 
(similar  to  the  BQR-7  stave)  having  a 39-inch  spacing  between  the  elements.  Broadening 
was  accomplished  by  delaying  the  outer  two  elements  equal  amounts  with  respect  to  the 
center  element,  thus  resulting  in  a vertical  time  delay  taper  which  is  an  even  function  of  ele- 
ment location.  This  results  in  a less  coherent  addition  of  signals  arriving  in  the  form  of  a 
plane  wave  along  the  vertical  boresight,  thereby  producing  a reduction  in  vertical  directivity. 

The  propagation  loss  characteristics  used  for  the  study  were  those  which  are  published 
in  USNUSL  Technical  Memo  2111-020-67.  The  ocean  model  is  the  Norwegian  Sea  (during 
summer)  having  a depth  of  2000  fathoms,  a layer  depth  of  100  feet,  and  sea  state  of  either 
2 or  4.  The  spectral  characteristics  of  the  target  model  were  also  taken  from  this  docu- 
ment. 
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2.  THEORETICAL  BEAM  PATTERN 

The  theoretical  beam  pattern,  as  implied  in  this  report,  is  the  normalized  plot  of  only 
that  portion  of  the  array  output  that  is  due  to  signal.  It  also  includes  the  effects  of  timing 
errors. 

If  the  instantaneous  output  of  the  array  is  considered  as  some  voltage,  the  average  value 
of  the  voltage  squared  is  proportional  to  the  average  power  out  of  the  array.  If  the  instan- 
taneous output  of  the  array  is  considered  as  some  digital  number,  the  average  value  of  the 
number  squared  is  the  variance  of  the  distribution  of  numbers  out  of  the  array.  In  either 
case,  there  will  be  some  power  (or  variance)  out  of  the  array  when  no  signal  i"  present  due 
to  noise.  The  signal  addition  will  cause  this  output  power  (or  variance)  to  change.  It  is  only 
this  change  in  output  which  is  of  significance  in  the  theoretical  beam  pattern. 

For  any  fixed  target  strength,  the  change  in  output  power  will  vary  with  the  angle  of 
arrival  of  the  plane  wave  from  the  target.  A three-dimensional  plot  of  how  this  change  in 
output  varies  with  the  angle  of  arrival,  and  normalized  to  a peak  of  unity,  defines  the  theo- 
retical beam  pattern. 

To  define  the  theoretical  beam  pattern  in  mathematical  terms,  let 
f.(t)  = total  output  of  element  i at  time  t 
T.  = applied  time  delay  in  output  of  element  i 
W.  = arbitrary  weighting  applied  to  output  of  element  i 
f^(t)  = total  output  of  array  at  time  t 
N = number  of  elements  used  in  array 

Since  the  array  output  is  formed  by  adding  the  delayed  outputs  of  each  element,  it  follows 

that 

N 

i=l 

The  output  of  each  element  is  composed  of  a noise  portion,  and  a signal  portion. 
f^.(t)  = noise  portion  of  output  of  element  i at  time  t 
f .(t)  = signal  portion  of  output  of  element  i at  time  t 
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then 

Vt)=fsi(t)+yt)  (5) 

Consider  some  hypothetical  omnidirectional  element  located  at  an  arbitrary  reference 
point.  Let 

f (t)  = output  of  hypothetical  reference  element  at  time  t due  to  a plane  wave  signal 
® from  a distant  point  source  target 

t.  --  transit  time  of  plane  wave  front  from  reference  point  to  element  i 

w.  = relative  response  of  element  i compared  to  hypothetical  standard;  Wj  can  be 
made  a function  of  angle  of  arrival  to  account  for  element  directivity 

From  the  above  definitions,  it  follows  that 

f .(t)  = w.  f (t  - t.)  (6) 

si' ' i s'  i'  ' ^ 

Squaring  Equation  (4),  taking  the  average  value,  and  combining  with  Equations  (5)  and 
(6)  leads  to  the  average  variance  (or  power)  out  of  the  array 


N N 
1=1  i=l 


A bar  over  a quantity  indicates  time  average  and  is  assumed  to  be  the  same  as  an 
ensemble  average. 

Since  the  signal  and  noise  are  independent  of  each  other,  the  average  value  of  cross 
product  terms  between  signal  and  noise  in  Equation  (7)  is  zero.  Equation  (7)  thus  simpli- 
fies to 


N N 


% - ‘i  - T()  - ‘j  - Tj) 

i=l  j=l 

+ W.  W.  f .(t  - T.)  f .(t  - T.) 
i j ni'  i'  nj'  ]' 


Notice  that  the  first  term  in  the  above  accounts  only  for  signal,  while  the  second  is 
only  noise.  In  keeping  with  the  definition  of  the  theoretical  beam  pattern  presented  earlier, 
interest  lies  only  in  the  signal  portion  of  this  output. 


a 

J 


[ 

r 

I 

L 

E 

[ 

[ 

r 

L 

E 


6 


CONFIDENTIAL 


CONFIDENTIAL 


The  time  delays  T.  are  chosen  such  that  when  the  signal  is  on  boresight,  T.  = -t..  On 

boresight  then,  the  time  averaged  quantity  reduces  to  f ^ (t) . Normalizing  the  signal  portion 

s 

of  Equation  (8)  to  this  peak  value  on  boresight  and  recognizing  the  normalized  autocorrela- 
tion function  of  the  signal  as  indicated  in  Equation  (9)  leads  to  Equation  (10), 


where 


P3(-ti-Ti-.yT.)  = 


(9) 


(10) 


p = theoretical  beam  pattern  of  linear  beamformer 


Finding  the  argument  of  the  autocorrelation  coefficient  in  terms  of  the  geometry  of  the 
array  still  remains.  Let 


r.  = vector  from  the  arbitrarily  chosen  reference  point  for  the  hypothetical  standard 
transducer  discussed  above  to  element  i 

c = magnitude  of  the  velocity  of  propagation  of  the  plane  wave  front  (signal) 

'c  = unit  vector  in  the  direction  of  propagation  of  the  plane  wave  front 


Using  conventional  notation  for  a vector  dot  product,  it  follows  that  the  transit  time 
from  the  reference  point  to  element  i is 


t.  = (-)  7. 

1 'c ' I 


(11) 


Let 

T - unit  vector  in  the  direction  chosen  for  boresight 

It  follows  from  the  definition  of  c and  ^ that  when  a target  is  on  boresight,  vectors  ^ 
and  T coincide.  Consequently,  the  exact  time  delay  required  for  beamforming  is  the  nega- 
tive of  the  transit  time  given  by  Equation  (II)  when  c is  replaced  by  s. 
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The  applied  time  delay  (T|)  may  deviate  slightly  from  this  ideal  value  due  to 
quantizing  errors.  The  error  in  time  delay  actually  used  will  then  be  given  by 


ei  = -tT|t(i)r:-T| 


The  argument  of  the  autocorrelation  coefficient  in  Equation  (10)  can  now  be  found  by 
combining  Equations  (11)  and  (12)  and  factoring.  Thus 

(-t.  - T.  + t.  + T.)  = (i)  (T.  - r.)  • (T-'H^)  + (e.  - e.)  (13) 

Notice  that  this  argument  is  independent  of  the  arbitrarily  chosen  origin.  It  depends 
only  on  the  vector  distance  between  elements  i and  j,  and  the  difference  between  boresight 
and  the  actual  angle  of  arrival. 

Equations  (10)  and  (13)  give  the  desired  theoretical  beam  pattern  by  letting  the  vector  'c 
scan  through  the  region  of  interest.  Arbitrary  spectra  are  accounted  for  by  appropriately 
specifying  the  normalized  autocorrelation  function  p . If  the  power  density  spectral  distri- 
bution  of  a signal  is  known,  its  autocorrelation  function  can  be  derived  by  a Fourier  trans- 
form relationship.  Element  directivity,  amplitude  weighting  to  shape  beams,  and  opaque 
arrays  can  be  accounted  for  by  the  factors  W.  w^.  (These  are  unity  for  omnidirectional 
elements  in  a transparent  array.)  Errors  in  timing  are  accounted  for  by  e.. 

When  computing  the  theoretical  beam  patterns,  it  has  been  assumed  that  the  element 
directivity  is  unity,  and  that  the  amplitude  weighting  is  uniform.  Timing  errors  which  are 
included  in  the  theoretical  beam  patterns  are  due  to  two  sources.  One  type  of  timing  error 
is  caused  by  time  delay  quantization.  The  second  type  of  error  results  from  using  a single 
comnromise  set  of  time  delays  for  two  or  more  classes  of  submarines. 


CONFIDENTIAL 


I V j 

t- 


CONFIDENTIAL 


SECTION  III 
RESULTS 

1.  VERTICAL  BEAM  COVERAGE  STUDY 
a.  Theoretical  Results 

One  of  the  objectives  of  this  study  was  to  determine  whether  the  detectability  of 
targets  could  be  improved  by  depressing  or  broadening  the  vertical  beam.  However,  the 
term  "detectability"  is  arbitrary  since  there  is  no  clear  cut  definition  of  what  is  required 
to  constitute  a target  detection.  The  approach  which  was  taken  here  was  to  calculate  the 
SNR  at  the  detector  output.  This  performance  measure  can  then,  at  least  in  theory,  be 
related  to  the  probability  of  detection.  Calculation  of  the  probability  of  detection,  however, 
is  not  considered  within  the  scope  of  this  report. 

The  computer  program  which  has  been  developed  calculates  the  SNR  for  a given 
input  target  range.  The  set  of  target  ranges  which  was  utilized  is  1.0,  1.6,  2.5,  4.0,  6.3, 
8.0,  10.0,  16,0,  25.0,  40.0,  and  50.0  kiloyards.  This  choice  of  ranges  was  used  since  the 
propagation  loss  data  for  these  parameters  has  been  published.  To  present  the  computer 
output  in  a useful  format  for  the  present  study,  the  SNR  was  calculated  for  these  ranges  at 
1*  increments  in  vertical  steering  angle  ranging  from  0 through  45®  depression  angle.  The 
data  was  then  linearly  interpolated  to  determine  the  contours  of  constant  SNR  on  a vertical 
depression  angle  versus  range  plot. 

Numerous  inputs  are  required  in  order  to  calculate  the  signal  power  spectrum  and 
the  noise  power  spectrum.  For  this  purpose  it  was  assumed  that  the  input  signal  spectrum 
was  weighted  by  the  present  BQR-7  pre-emphasis  filter  and  that  the  power  response  char- 
acteristic of  the  transducer  was  uniform.  Furthermore,  the  own  ship's  self -noise  level  was 
modeled  after  the  SSB(N)  616  platform  moving  at  a speed  of  five  knots. 

All  SNR  contour  plots  shown  in  this  section  of  the  report  are  evaluated  on  the  basis 
of  an  output  filter  integration  time  of  0.5  second.  In  other  words,  the  term  in  the  de- 
nominator of  Equation  (1)  has  been  set  equal  to  unity.  If  a filter  integration  time  other 
than  0.5  second  is  used,  the  SNR  is  theoretically  shifted  by  10  log^^  (2T).  If  an  integi'ation 
time  of  12.8  seconds  is  used,  as  is  anticipated  for  the  BQR-7  DIMUS,  the  SNR  curve  is 
shifted  by  about  14  dB. 
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Three  target  models  for  which  the  spectral  characteristics  were  available  were 
used  to  calculate  the  signal  power  spectrum.  These  were  a 10-knot  snorkler,  a Pkipjack 
[SS(N)  585] , and  a Permit  [SS(N)  594]  submarine.  Figures  1 through  3 are  the  respective 
SNR  contour  plots  for  these  three  targets  in  the  absence  of  vertical  beam  broadening.  The 
results  for  these  plots  are  based  on  the  assumption  that  the  target  and  the  receiver  are  lo- 
cated on  the  same  side  of  the  layer.  The  plots  are  repeated  in  Figures  4 through  6 for  the 
case  where  the  target  and  the  receiver  are  located  on  opposite  sides  of  the  layer.  Sea  state 
2 is  assumed  in  each  case. 

As  expected,  the  high  SNR  contours  decrease  in  range  as  the  vertical  beam  is 
depressed.  This  occurs  since  the  energy  in  the  surface  duct  (which  is  the  primary  propa- 
gation mode  at  short  range)  enters  the  system  off  the  main  response  axis  of  the  vertical 
beam  pattern.  For  noisy  targets  such  as  the  10 -knot  snorkler,  a definite  range  improvement 
is  noted  by  depressing  the  vertical  beam  by  9 or  10°.  It  should  be  emphasized  that  this  im- 
provement is  totally  due  to  the  fact  that  the  target  is  extremely  noisy  and  that  the  signal 
energy  in  the  bottom  bounce  propagation  mode  is  sufficient  to  yield  a high  SNR  at  long  ranges. 
The  results  indicate  that  the  24-dB  contour  extends  beyond  the  50-kyd  range  for  an  integra- 
tion time  of  0.5  second  regardless  of  whether  the  target  and  receiver  are  located  on  the 
same  or  opposite  sides  of  the  layer.  With  an  integration  time  of  12.8  seconds  each  of  the 
contours  shown  is  theoretically  improved  by  an  additional  14  dB. 

For  quiet  targets  such  as  the  Permit  submarine,  an  improvement  in  range  is  not 
so  obvious.  It  tends  to  be  limited  by  the  assumptions  with  respect  to  the  environment, 
directivity,  and  the  improvement  resulting  from  an  increased  integration  time.  If  there  is 
no  loss  in  the  vertical  directivity  due  to  depressing  the  beam  and  the  full  14  dB  of  SNR  im- 
provement is  realized  by  the  BQR-7  DIMUS  12.8  seconds  integration  time,  the  0-dB  contour 
would  fall  between  the  present  -12  and  -18  dB  contour  in  Figure  3.  This  would  then  result 
in  improved  range  performance  for  the  0-dB  contour  if  the  beam  were  depressed  10  to  15°. 

It  should  be  emphasized  however,  that  if  in  practice  the  SNR  improvement  is  closer  to  6 dB, 
then  the  0-dB  contour  gives  considerably  worse  range  performance  than  it  would  if  the 
vertical  beam  had  not  been  depressed  at  all.  The  same  general  conclusions  can  be  made 
with  respect  to  the  case  where  the  target  and  the  receiver  are  situated  on  opposite  sides  of 
the  layer  (see  Figure  6).  If  the  theoretical  improvement  due  to  time  averaging  12.8  seconds 
could  be  achieved,  the  best  0-dB  range  contour  would  Ix'  obtained  by  depressing  the  beam 
between  20  and  25°. 
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SNR  Contours  for  a Permit  Submarine  with  Target  and  Receiver 
on  Same  Side  of  flayer 
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Figure  5.  SNR  Contours  for  a Skipjack  Submarine  with  Target  and  Receiver 
on  Opposite  Sides  of  Layer 


CONFIDENTIAL 


CONFIDENTIAL 

SNR  contour  plots  including  the  effects  of  vertical  beam  broadening  using  the  Permit 
submarine  for  target  radiated  noise  were  computed.  This  target  model  was  chosen  since  it 
has  the  least  favorable  (quietest)  characteristics  for  detection.  Before  presenting  the 
results,  a few  comments  are  in  order  to  clarify  the  method  used  to  achieve  the  vertical 
beam  broadening. 

Computerized  implementation  of  vertical  beam  broadening  was  accomplished  by 
appropriate  adjustment  of  the  vertical  beam  shape  factor  of  Equation  (2).  The  vertical  beam 
pattern  was  assumed  to  be  independent  of  the  horizontal  steering  angle.  Furthermore,  the  ver- 
tical pattern  was  assumed  to  be  equivalent  to  that  of  a three-element  line  array.  Broadening 
was  achieved  by  delaying  the  outer  two  elements  by  equal  amounts.  Originally  it  was  in-  | 

tended  to  accomplish  beam  broadening  by  computing  the  applied  time  delays  on  the  basis  of  * 

an  erroneous  velocity  of  sound  which  is  higher  than  that  experienced  at  sea.  Since  it  was 
more  feasible  however,  from  the  standpoint  of  existing  computer  programs  to  implement 
the  broadening  procedure  as  outlined  above,  the  former  approach  was  selected.  This  method 

I 

also  has  the  flexibility  of  varying  the  depression  angle  and  the  broadening  independently. 

To  illustrate  that  vertical  beam  broadening  can  be  achieved  ty  use  of  this  method, 
the  vertical  beam  patterns  for  the  BQR-7  array  were  calculated  using  various  amounts  of 
delay  for  the  upper  and  lower  layers  (or  rows)  of  elements.  The  results  are  shown  in 
Figures  7 through  13.  Each  of  these  figures  shows  the  vertical  beam  pattern  for  seven 
specific  beams.  Beams  1 through  3 are  directed  toward  the  stern  of  the  ship,  Beam  12  is 
directed  toward  the  port  side,  and  Beams  29  through  31  are  directed  forward.  Note  that  a 
separate  abscissa  is  drawn  for  each  beam.  Each  beam  is  normalized  so  that  the  0-dB  point  | 
represents  completely  coherent  addition  of  the  input  signal  (a  theoretical  maximum) . Any  : 

loss  is  therefore  caused  by  the  time  delay  errors  required  for  broadening  or  time  delay 
errors  due  to  quantization.  Ticks  are  drawn  at  3-dB  intervals  from  the  peak  response  of 
each  beam. 

Figure  7 represents  the  vertical  patterns  without  broadening.  Figures  8 through  13 
indicate  the  results  obtained  by  applying  delays  up  to  375  jis  in  steps  of  62.5  ns.  Note  spe- 
cifically how'  the  beam  is  broadened  by  increasing  amounts  of  delay.  Further  note  that 
broadening  is  accompanied  by  a noticeable  reduction  in  gain  along  boresight.  I 

j 
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Vertical  Beam  Patterns  Without  Broadening 
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Figure  8.  Vertical  Beam  Patterns  Obtained  by  Delaying  Upper  and  Lower  Rows 
of  Elements  by  62.5  ns 
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Figure  10.  Vertical  Beam  Patterns  Obtained  by  Delaying  Upper  and  Lower  Rows 
of  Elements  by  187.5  ns 
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Figures  14  through  16  are  the  SNR  contours  for  various  amounts  of  broadening  with 
the  target  and  receiver  on  the  same  side  of  the  layer.  Figure  14  indicates  the  results  ob- 
tained by  applying  a 62.5-jns  delay  to  each  of  the  outer  elements  of  the  stave  while  Figures  15 
and  16  are  the  results  for  appl3dng  125-  and  187.5-/is  delays,  respectively.  Figures  17 
through  19  indicate  the  results  for  identical  amounts  of  broadening  but  with  the  target  and 
receiver  on  opposite  sides  of  the  layer. 

The  particular  delays  chosen  to  implement  vertical  beam  broadening  for  the  SNR 
contours  were  selected  to  give  a specified  amount  of  phase  shift  between  the  outer  and  center 
elements.  The  resulting  phase  shifts  at  the  band  center  frequency  are  7,  and  ^ radians. 

O t:  O 

One  observes  the  same  trends  for  broadening  as  were  observed  for  depression  only. 
Low  SNR  contours  show  improved  range  performance.  The  range  improvement  for  low  SNR 
contours  is  achieved  however  at  the  expense  of  reduced  performance  at  higher  SNR  contours. 
The  detectability  of  quiet  targets,  such  as  the  Permit,  at  shorter  ranges  has  been  impaired. 

Figures  20  and  21  show  the  results  for  a Permit  target,  no  vertical  beam  broaden- 
ing and  sea  state  4,  Note  from  these  curves  that  even  if  a 12-dB  Improvement  were  achieved 
through  Integration,  beam  depression  would  result  in  lower  ranges  for  the  0-dB  SNR  contour. 
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b.  Results  from  Sea  Test  Data 

An  attempt  was  made  to  verify  the  theoretical  results  by  beamforming  taped  data 
recorded  at  sea.  The  tape  is  a recording  of  the  BQR-7  preamplifier  outputs  aboard  the 
SSB(N)  657,  Francis  Scott  Key.  The  USS  Angler  served  as  the  target  ship.  Although  Key 
was  on  her  first  exercise  following  commissioning,  the  performance  of  her  sonar  appeared 
to  be  excellent. 

Tape  No.  10  was  utilized  for  the  present  experimental  verification,  and  taken  while 
Angler  circled  Key  at  a range  of  about  2000  yards.  Angler  was  snorkeling  while  Key, 
situated  at  periscope  depth,  maintained  a 3-knot  straight  course.  Sea  states  2 and  3 existed. 
The  water  depth  was  approximately  1650  fathoms.  The  velocity  of  sound  was  estimated  at 
4975  feet  per  second. 

Figures  22  through  24  illustrate  the  correlation  patterns  obtained  by  hard-limiting, 
sanipling,  and  beamforming  the  hydrophone  outputs.  The  ordinate  values  on  each  of  the  plots 
correspond  to  a digital  number  proportional  to  the  beamformer  output. 

Figure  22  illustrates  plots  of  all  62  of  the  BQR-7  DIMUS  beams.  The  velocity  of 
sound  used  to  compute  the  delays  for  beamforming  this  data  is  5120  feet  per  second,  which 
is  145  feet  per  second  higher  than  estimated  actual  velocity.  The  steering  angles  used  to 
beamform  the  data  were  chosen  on  the  basis  of  giving  the  best  coverage  for  the  Key  hull  con- 
figuration. 

Figures  23  and  24  are  plots  of  the  outputs  of  four  of  the  aft  beams  (Beams  2 through 
5)  and  four  forward  beams  (Beams  28  through  31)  using  a sound  velocity  of  4950  feet  per 
second.  Performance  of  only  those  beams  most  sensitive  to  the  velocity  of  sound  were 

I 

verified.  This  was  primarily  done  to  conserve  computer  time  since  the  beamforming  runs 
are  quite  lengthy. 

If  the  theoretical  results  are  to  be  verified  by  actual  sea  test  data,  then  the  response 
curves  of  Figures  23  and  24  should  show  better  performance  than  the  corresponding  curves 
illustrated  in  Figure  22.  This  should  be  true  especially  since  the  target  was  at  close  range 
where  the  surface  duct  is  the  primary  propagation  mode.  Consequently,  beam  depression 
or  vertical  beam  broadening  should  hinder  rather  than  enhance  detection  of  the  target.  Upon 
comparing  the  response  curves  of  the  respective  beams,  it  is  observed  that  the  best  perform- 
ance is  obtained  by  calculating  time  delays  using  a sound  velocity  of  4950  feet  per  second. 
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2.  COMPROMISE  BEAMFORMER  FEASIBILITY  STUDY 

It  was  the  Liteiit  of  this  phase  of  the  study  to  determine  whether  a suitable  compromise 
beamformer,  or  perhaps  two  compromise  beamformers,  could  be  designed  to  accommodate 
the  following  four  classes  of  submarines. 

SS(N)  594* 

SS(N)  671 
SSB(N)  598 
SSB(N)  608** 

The  definition  of  a compromise  beamformer,  as  implied  here,  is  that  it  must  provide 
suitable  fixed  steering  delays  to  perform  the  beamforming  function  for  more  than  one  class 
of  submarines.  The  selection  of  the  time  delays  for  the  compromise  beamformer  must 
therefore  be  made  on  the  basis  of  optimizing  the  combined  performance  of  more  than  one 
class  of  submarines. 

To  provide  a standard  of  comparison,  the  theoretical  horizontal  beam  patterns  were 
calculated  for  each  of  the  four  classes  of  submarines.  These  are  the  patterns  which  are 
obtained  by  computing  the  steering  delays  for  each  class  of  submarines  assuming  a sound 
velocity  of  4950  feet  per  second.  These  delays  are  then  quantized  to  the  nearest  62-/iS 
interval.  The  pattern  is  subsequently  generated  using  these  delays  in  conjunction  with 
Equation  (10).  For  the  purpose  of  generating  the  theoretical  patterns,  a flat  spectrum 
having  a bandwidth  of  1000  Hz  and  centered  at  1000  Hz  was  assumed.  This  results  in  a 
normalized  autocorrelation  function  as  depicted  in  Figure  25. 

The  theoretical  beam  patterns  are  plotted  in  Figures  26  through  29.  Only  31  of  the  82 
beams  are  plotted  since  the  other  31  beams  are  identical  to  those  shown  (because  of  the 
symmetry  of  the  array).  Note  that  in  each  of  the  theoretical  patterns  three  of  the  beams 
are  appreciably  degraded.  This  results  from  the  fact  that  some  of  the  steering  delays  for 
these  beams  are  larger  than  the  available  delays  of  the  beamformer.  Consequently,  the 
output  from  these  elements  is  added  with  a delay  which  is  less  than  required.  The  delays 
however,  are  selected  so  that  there  will  be  positive  correlation  with  the  output  of  the  re- 
maining elements.  This  is  achieved  by  selecting  a proper  time  delay  based  on  the 


* 'I'he  results  for  the  SS(N)  594  are  also  applicable  to  the  SS(N)  637. 

**  The  results  for  the  SSB(N)  608  are  also  applicable  to  the  SSB(N)  616  and  the  SSB(N)  640. 
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Figure  25.  Normalized  Autocorrelation  Function  for  Assumed  Signal  Spectrum 
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Figure  26.  Theoretical  horizontal  Beam  Patterns  for  SS(N)  594 
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Figure  28,  Theoretical  Horizontal  Beam  Patterns  for  SSB(N)  608 
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normalized  autocorrelation  function.  In  practice  the  elements  which  are  improperly  de- 
layed degrade  the  beam  output  very  little  since  the  elements  which  are  affected  have 
very  little  response  in  the  boresight  direction  of  the  beams  involved.  The  peak  responses 
for  the  theoretical  beam  patterns  for  the  four  classes  of  submarines  are  tabulated  in  Table  1. 


TABLE  1 

PEAK  RESPONSE  OF  THEORETICAL  BEAM  PATTERNS 


Gain  (dB) 

SS(N)  594 

SS3(N)  598 

SSB(N)  608 

SS(N)  671 

Beam 

Number 

1 

-1.69 

-1.51 

-1.69 

-1.78 

2 

-1.04 

-1.04 

-1.05 

-1.17 

3 

-0.56 

-0.56 

-0.56 

-0.56 

4 

-0.03 

-0.05 

-0.12 

-0.15 

5 

-0.04 

-0.04 

-0.04 

-0.04 

6 

-0,04 

-0.04 

-0.04 

-0.05 

7 

-0.03 

-0.04 

-0.04 

-0.04 

8 

-0.04 

-0.04 

-0.04 

-0.04 

9 

-0.04 

-0.04 

-0.04 

-0.05 

10 

-0.04 

-0.04 

-0,04 

-0.03 

11 

-0.04 

-0.04 

-0.03 

-0.04 

12 

-0.04 

-0.04 

-0.03 

-0.04 

13 

-0.04 

-0.04 

-0.04 

-0.04 

14 

-0.04 

-0.03 

-0.04 

-0.03 

15 

-0.04 

-0.03 

-0.03 

-0.03 

16 

-0.04 

-0.04 

-0.04 

-0.04 

17 

-0.05 

-0.04 

-0.04 

-0.03 

18 

-0.03 

-0.04 

-0.04 

-0.03 

19 

-0.04 

-0.05 

-0.04 

-0.05 

20 

-0.04 

-0.04 

-0.03 

-0.04 

21 

-0.03 

-0.04 

-0.04 

-0.04 

22 

-0.04 

-0.03 

-0.04 

-0.04 

23 

-0.04 

-0.04 

-0.04 

-0.04 

24 

-0.04 

-C.04 

-0.04 

-0.05 

25 

-0.05 

-0.04 

-0.04 

-0.04 

26 

-0.04 

-0.04 

-0.04 

-0.04 

27 

-0.04 

-0.05 

-0.03 

-0.04 

28 

-0.04 

-0.04 

-0.04 

-0.03 

29 

-0.05 

-0.04 

-0.04 

-0.05 

30 

-0.04 

-0.04 

-0.04 

-0.04 

31 

-0.04 

-0.04 

-0.03 

-0.04 
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Compromise  beamformers  were  investigated  for  various  combinations  of  array  config- 
urations (classes  of  submarines).  The  delay  values  for  the  compromise  beamformer  were 
calculated  on  the  basis  of  the  following  method.  The  average  delay  value  per  element 
was  computed  for  each  beam  and  class  of  submarine.  The  average  delay  value  was  then 
adjusted  so  that  the  corresponding  beam  of  each  class  of  submarine  possessed  the  same 
average  delay  per  element.  The  adjustment  was  made  by  adding  or  subtracting  a constant 
to  each  delay  value  for  each  beam  (different  for  each  beam  but  the  same  for  each  element). 
This  does  not  alter  the  pattern  since  it  is  only  the  delay  difference  between  elements  which 
is  instrumental  in  producing  the  pattern.  Thereafter,  the  mean  value  of  the  delays  of  the 
corresponding  elements  in  corresponding  beams  was  computed  to  determine  the  compromise 
delay  values.  Any  constant  delay  which  had  been  added  or  subtracted  to  adjust  the  mean 
delay  per  element  was  then  reinstated. 

Figures  30  through  33  illustrate  the  beam  patterns  for  the  four  classes  of  submarines 
if  a single  compromise  beamformer  was  designed  to  handle  all  four  classes.  Table  2 lists 
the  corresponding  peak  responses  for  each  class  submarine  using  the  four-array  compro- 
mise and  plots  of  these  losses  are  shown  in  Figure  34. 

Figures  35  through  37  are  the  beam  patterns  for  the  SS(N)  594,  SSB(N)  598,  and  the 
SSB(N)  608  using  a compromise  beamformer  for  these  three  classes  of  submarines. 

The  corresponding  losses  in  peak  response  are  tabulated  in  Table  3 and  a plot  of  these 
losses  is  shown  in  Figure  38.  Various  other  two-array  compromises  were  attempted.  The 
corresponding  beam  patterns  are  shown  in  Figures  39  through  51.  The  resulting  losses  in 
peak  response  are  tabulated  in  Table  4 and  plotted  in  Figures  52  through  54. 

The  "worst  case"  losses  in  peak  response  resulting  from  all  combinations  of  compro- 
mise delay  values  were  tabulated  on  a beam  to  beam  basis.  These  are  plotted  in  Figure  55. 

The  arrival  times  and  the  sampling  times  for  numerous  configurations  have  been  com- 
pilea  and  are  shown  in  Figures  56  through  63.  The  arrival  time  is  computed  by  subtracting 
the  delay  time  from  a positive  constant  greater  than  the  largest  required  delay.  The  re- 
lationship between  the  delay  times  and  the  sampling  times,  however,  is  not  so  obvious. 

Using  the  beamformer  design  as  given  in  "Final  Report,  AN/BQR-7,  All  Digital  Px’eformed 
Beam  Receiver  Study  (U)"  dated  September  30,  1965,  the  sampling  times  can  be  computed 
from  the  arrival  time  as  outlined  in  this  document. 
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Figure  30.  Horizontal  Beam  Patterns  of  SS(N)  594  Using  a Four-Array 
Compromise  Beamformer  Based  on  SS(N)  594,  SSB(N)  598, 
SSB(N)  608,  and  SS(N)67l  Configurations 
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Figure  31.  Horizontal  Beam  Patterns  of  SSB(N)  598  Using  a Four-Array 
Compromise  Beamformer  Based  on  SS(N)  594,  SSB(N)  598, 
SSB(N)608,  and  SS(N)67l  Configurations 
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Figure  32.  Horizontal  Beam  Patterns  of  SSB(N)  608  Using  a Four-Array 
Compromise  Beamformer  Based  on  SS(N)  594,  SSB(N)  598, 
SSB(N)  608,  and  SS(N)  671  Configurations 
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Figure  33.  Horizontal  Beam  Patterns  of  SS(N)  671  Using  a Four-Array 
Compromise  Beamformer  Based  on  SS(N)  594,  SSB(N)  598, 
SSB(N)  608,  and  SS(N)  671  Configurations 
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TABLE  2 

PEAK  RESPONSE  VALUES  FOR  FOUR-ARRAY  COMPROMISE 


Gain  (dB) 

1 

Oi 

lO 

§ 

CO 

to 

SSB{N)  598 

SSB(N)  608 

SS(N)  671 

Beam 

Number 

1 

-1.73 

-1.83 

-1.71 

-1.85 

2 

-1.09 

-1.28 

-1.08 

-1.22 

3 

-0.60 

-0.82 

-0.60 

-0.77 

4 

-0.13 

-0.25 

-0.14 

-0.24 

5 

-0.15 

-0.32 

-0.11 

-0.28 

6 

-0.12 

-0.27 

-0.10 

-0.37 

7 

-0.14 

-0.34 

-0.08 

-0.42 

8 

-0.11 

-0.32 

-0.02 

-0.46 

9 

-0.10 

-0.41 

-0.04 

-0.57 

10 

-0.07 

-0.39 

-0.07 

-0,62 

11 

-0.12 

-0.50 

-0.06 

-0.57 

12 

-0.07 

-0.58 

-0.08 

-0.58 

13 

-0.12 

-0.54 

-0.12 

-0.49 

14 

-0,06 

-0.49 

-0.08 

-0.64 

15 

-0.10 

-0.51 

-0.10 

-C.69 

16 

-0.13 

-0.62 

-0.08 

-0.60 

17 

-0.08 

-0.53 

-0.14 

-0.66 

18 

-0.08 

-0.59 

-0.14 

-0.72 

19 

-0.10 

-0.54 

-0.13 

-0.57 

20 

-0.11 

-0.64 

-0.10 

-0.52 

21 

-0.11 

-0.57 

-0.11 

-0.42 

22 

-0.10 

-0.52 

-0.09 

-0.26 

23 

-0.09 

-0.44 

-0.17 

-0.28 

24 

-0.16 

-0.40 

-0.16 

-0.29 

25 

-0.07 

-0.44 

-0.16 

-0.26 

26 

-0.14 

-0.38 

-0.21 

-0.19 

27 

-0.15 

-0.45 

-0.18 

-0.17 

28 

-0.11 

-0.54 

-0.20 

-0.12 

29 

-0.09 

-0.66 

-0.23 

-0.12 

30 

-0.20 

-0.67 

-0.24 

-0.20 

31 

-0.12 

-0.72 

-0.31 

-0.21 
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Figure  34.  Plot  of  Peak  Response  Values  for  Four-Array  Compromise 
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Figure  35.  Horizontal  Beam  Patterns  of  SS(N)  594  Using  a Three-Array 
Compromise  Beamformer  Based  on  SS(N)  594,  SSB(N)  598, 
and  SSB(N)  608  Configurations 
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Figure  36.  Horizontal  Beam  Patterns  of  SSB(N)  598  Using  a Three-Array 
Compromise  Beamformer  Based  on  SS(N)  594,  SSB(N)  598, 
and  SSB(N)  608  Configurations 
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TABLE  3 

PEAK  RESPONSE  VALUES  FOR  THREE-ARRAY  COMPROMISE 


Gain  (dB) 

SS(N)  594 

SSB(N)  598 

SSB(N)  608 

Beam 

Number 

1 

-1.70 

-1.82 

-1.67 

2 

-1.09 

-1.26 

-1.09 

3 

-0.65 

-0.69 

-0.62 

4 

-0.13 

-0.22 

-0.08 

5 

-0.10 

-0.26 

-0.09 

6 

-0.13 

-0.19 

-0.10 

7 

-0.12 

-0.23 

-0.06 

8 

-0.13 

-0.16 

-0.12 

9 

-0.12 

-0.17 

-0.12 

10 

-0.10 

-0.20 

-0.13 

11 

-0.13 

-0.16 

-0.16 

12 

-0.08 

-0.27 

-0.10 

13 

-0.09 

-0.22 

-0.17 

14 

-0.09 

-0.22 

-0.17 

15 

-0.09 

-0.25 

-0.16 

16 

-0.09 

-0.30 

-0.14 

17 

-0.09 

-0.32 

-0.17 

18 

-0.10 

-0.33 

-0.19 

19 

-0.13 

-0.36 

-0.16 

20 

-0.10 

-0.43 

-0.13 

21 

-0.14 

-0.38 

-0.15 

22 

-0.13 

-0.31 

-0.24 

23 

-0.15 

-0.35 

-0.20 

24 

-0.17 

-0.41 

-0.19 

25 

-0.13 

-0.38 

-0.17 

26 

-0.16 

-0.44 

-0.19 

27 

-0.18 

-0.40 

-0.24 

28 

-0.14 

-0.52 

-0.17 

29 

-0.24 

-0.37 

-0.34 

30 

-0.22 

-0.55 

-0.30 

31 

- — — 

-0.19 

-0.51 

-0.39 
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Figure  38,  Plot  of  Peak  Response  Values  for  Three-Array  Compromise 
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Figure  39.  Hor'-sontal  Beam  Patterns  of  SS(N)  594  Using  a Two-Array  Compromise 
Beamformer  Based  on  SS(N)  594  and  SSB(N)  598  Configurations 
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Figure  40.  Horizontal  Beam  Patterns  of  SSB(N)  598  Using  a Two-Array  Compromise 
Beamformer  Based  on  SS(N)  594  and  SSB(N)  598  Configurations 
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Figure  41.  Horizontal  Beam  Patterns  of  SSB(N)  608  Using  a Two-Array  Compromise 
Beamformer  Based  on  SSB(N)  608  and  SS(N)  671  Configurations 
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Figure  43.  Horizontal  Beam  Patterns  of  SS(N)  594  Using  a Two-Array  Compromise  I 

Beamformer  Based  on  SS(N)  594  and  SSB(N)  608  Configurations  | 
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Figure  44,  Horizontal  Beam  Patterns  of  SSB(N)  608  Using  a Two-Array  Compromise 
Beamformer  Based  on  SS(N)  594  and  SSB(N)  608  Configurations 
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Figure  45.  Horizontal  Beam  Patterns  of  SSB(N)  598  Using  a Two-Array  Compromise 
Beamformer  Based  on  SSB(N)  598  and  PS(N)  671  Configurations 
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Figure  46.  Horizontal  Beam  Patterns  of  SS(N)  671  Using  a Two-Array  Compromise 
Beamformer  Based  on  SSB(N)  598  and  SS(N)  671  Configurations 
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Figure  47.  Horizontal  Beam  Patterns  of  SS(N)  594  Using  a Two-Array  Compromise 
Beamformer  Based  on  ?S(N)  594  and  SS(N)  671  Configuration 
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Figure  48.  Horizontal  Beam  Patterns  of  SS(N)  671  Using  a Two-Array  Compromise 
Beamformer  Based  on  SS(N)  594  and  SS(N)  671  Configurations 
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Figure  49.  Horizontal  Beam  Patterns  of  SSB(N)  608  Using  a Two-Array  Compromise 
Beamformer  Based  on  SSB(N)  598  and  SSB(N)  608  Configurations 
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Figure  50.  Horizontal  Beam  Patterns  of  SSB(N)  598  Using  a Two-Array  Compromise 
Beamformer  Based  on  SSB(N)  598  and  SSB(N)  608  Configurations 
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Figure  52.  Plot  of  Peak  Response  Values  for  a Two-Array  Compromise 
Between  the  SS(N)  594  and  SS(N)  671  and  a Two-Array  Com- 
promise Between  the  SSB(N)  608  and  SSB(N)  598  Configurations 
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Figure  53,  Plot  of  Peak  Response  Values  for  a Two-Array  Compromise 
Between  the  SS(N)  594  and  SSB(N)  608  and  a Two-Array  Com- 
promise Between  the  SSB(N)  598  and  SS(N)  671  Configurations 


CONFIDENTIAL 


ESPONSE  (dB) 


CONFIDENTIAL 


o SS(N)  594 
□ SSB(N)59B 
A SSB(N)608 
V SS(N)671 


BEAM  STEERING  ANGLE  (DEGREES) 


Figure  51.  Plot  of  Peak  Response  Values  for  a Two-Array  Compromise 
Between  the  SSB(N)  608  and  SS(N)  671  and  a Two-Array  Com- 
promise Between  the  SS(N)  594  and  SSB(N)  598  Configurations 
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Figure  55,  Worst  Case  Loss  in  Peak  Response  Resulting  from  Various  Combinations 
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Figure  56.  Arrival  Times  for  SS(N)  594  Class 
Submarine  (Sheet  1) 
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SECTION  IV 
CONCLUSIONS 

1.  VERTICAL  BEAM  COVERAGE  STUDY 

It  is  difficult  to  make  a clear-cut  decision  as  to  whether  vertical  beam  broadening  and 
depression  are  beneficial  to  the  overall  detection  problem.  For  noisy  targets,  such  as  a 
snorkler,  the  detection  range  is  improved  somewhat  by  beam  depression.  For  quiet  targets, 
such  as  the  Permit  submarine,  the  detection  range  may  be  impaired.  Beam  depression  and 
vertical  beam  broadening  is  effective  only  in  cases  where  array  gain  for  surface  duct  re- 
ception can  be  sacrificed  at  the  expense  of  improved  bottom  bounce  reception.  This  trade- 
off can  be  costly  in  terms  of  detectability  of  quieter  targets. 

It  is  evident  that  any  conclusions  are  very  sensitive  to  the  characteristics  of  the  model. 
Whereas  one  set  of  specified  inputs  might  dictate  a certain  degree  of  beam  depression  and 
broadening,  this  might  be  detrimental  to  target  detectability  for  another  set  of  Inputs,  To 
assure  the  detectability  of  quiet  targets  with  the  array  gain  provided  by  the  BQR-7  DIMUS, 
it  appears  that  vertical  beam  broadening  and  depression  is  not  desirable. 

2.  COMPROMISE  BEAMFORMER  FEASIBILITY  STUDY 

Several  of  the  classes  of  submarines  are  quite  compatible  for  the  purpose  of  designing 
a suitable  compromise  beamformer.  Chief  candidates  for  a suitable  two-array  compromise 
are  the  following  pairs. 

SS(N)  594.  and  SSB(N)  608 
SSB(N)  608  and  SS(N)  671 

Degradation  of  the  peak  response  for  most  beams  is  less  than  0,2  dB  utilizing  these 
pairings.  The  least  desirable  combinations  for  a two-array  compromise  are  the  following 
pairs. 

SSB(N)  598  and  SS(N)  671 
SSB(N)  598  and  SSB(N)  608 

The  degradation  for  these  pairings  is  considered  unnecessarily  large  when  compared 
to  other  two-array  compromises. 

The  largest  degradation  for  the  three-array  and  four-array  compromises  which  were 
considered  in  Section  III  are  approximately  0.5  and  0.8  dB,  respectively. 

nu 
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If  a peak  loss  of  0.8  dB  due  to  the  effects  of  compromising  the  delays  is  considered 
acceptable,  the  four-array  compromise  is  suitable.  If  a peak  loss  in  the  vicinity  of  0.5  dB 
is  acceptable,  a three-array  compromise  between  the  SS(N)  594,  SSB(N)  598,  and  the 
SSB(N)  608  is  recommended.  If  both  of  these  degradations  are  excessive,  the  best  pair  of 
two-array  compromises  which  include  all  four  classes  of  submarines  is  the  following. 

SSB(N)  608  and  SS(N)  671 
SS(N)  594  and  SSB(N)  598 

Tne  largest  degradation  for  these  pairings  is  approximately  0.3  dB. 


